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Abstract. We discuss matching control laws for undcractuated systems. We previously showed 
that this class of matching control laws is completely charactarized by a linear system of first order 
partial differential equations for one set of variables ( A) followed by a linear system of first order 
PDEs for the second set of variables (g, V). Here we derive a new first order system of partial 
differential equations that encodes all compatibility conditions for the A-equations. We give four 
examples illustrating different features of matching control laws. The last example is a system with 
two unactuated degrees of freedom that admits only basic solutions to the matching equations. There 
are systems with many matching control laws where only basic solutions are potentially useful. We 
introduce a rank condition indicating when this is likely to be the case. 
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1. Introduction . Effective procedures for designing control laws are very im- 
portant in nonlinear control theory. Explicit analytic formulae for control laws play a 
role similar to explicit solutions to differential equations. Such formulae exist in only 
a few special cases, but those that exist serve as simple models to develop and test 
more general techniques. 

In this paper we discuss a class of full state feedback control laws for underactu- 
ated systems. In Q we showed that this class of matching control laws is completely 
charactarized by a linear system of first order partial differential equations for one 
set of variables ( A) followed by a linear system of first order PDEs for the second 
set of variables {g, V). These equations always have a simple family of solutions 
which we call basic solutions. The system of equations for the first set of variables 
( A-equations) is overdetermined. Here we derive a new first order system of partial 
differential equations that encodes all compatibility conditions for the A-equations 
(we call these the i/-equations). If only one degree of freedom is unactuated, the solu- 
tions to all these systems of PDEs can be completely analyzed. It is often possible to 
get explicit formulae for the solutions to these equations. We also provide an example 
of a system with two unactuated degrees of freedom that has only basic solutions. 
There are systems with many matching control laws where only basic solutions are 
potentially useful. We write down a rank condition indicating when this is likely to 
be the case. 

During the last few years several researchers have investigated control laws in 
which the closed loop system assumes a certain structure. Numerous papers have 
been written on this subject, see [1] - [13] and the references therein. The control 



laws that form the subject of this present paper are described by equations (2.4) and 
(2.6). These equations were independently derived in | [lO| and ||^. The A-equations 



were first introduced in paper Even though the initial matching equations of 
[ p^ and Q form a highly nonlinear system of PDEs, introduction of the A variables 
triangulates the system. The system is triangulated in the sense that all solutions are 
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obtained by first solving first order linear equations for A and then solving first order 
linear equations for the remaining variables. 

This paper is organized as follows. Section 2 reviews matching control laws and 
the A-equations, and introduces the z^-equations. Section 3 specializes to systems with 
one unactuated degree of freedom. Sections 4, 5 and 6 contain examples illustrating 
three different features of matching control laws. The rank condition appears at the 
end of Section 5. Later we apply it in Section 6. In Section 7 we describe the final 
example of a system with two unactuated degrees of freedom. We show that this 
system has only basic matching control laws. 

2. Matching equations. We use the following notation. 

• n is the number of the degrees of freedom of the mechanical system 

• X = (x^ , . . . , x") are configuration variables denoting the position of the system, 

and X = (i^, . . . ,x") are the corresponding velocities 

• 9ij{^) is the mass matrix 

• V{x) is the potential energy 

• Ci{x,x) are the dissipation terms 

• Ui{x,x) are the control inputs 

Let m < n be the number of unactuated degrees of freedom. We will assume that 
degrees of freedom numbered 1 through m are unactuated and use indices a, 6, ... 
to indicate unactuated degrees of freedom. The indices i, j, ... will run from 1 to n. 
We adopt the convention of summation over the repeated indices . 
Given this, the equations of motion of the system are 

■ fc 9V 

grjX^ + [j k, r]x^ X + + = Wr , r = l,...,n, (2.1) 
where [ij, k] is the Christoffel symbol of the first kind. 

Our assumtion that the first m degrees of freedom are not actuated means that 

ui = ■■■ ^ it,,^ = 0. (2.3) 

We are looking for control laws Ui such that the closed loop system can be written in 
the form 

■ I - dV 

9rjx' + [j k, r]x^ X + Cr + =0, r=l,...,n. 



where [ij, k\ is defined as in ( |2.2| ) with g in place of g. Such a control law will be 
given by 

= ( [j k, e] - gu T [Jfc^] ) + {Ci- gu ^ C, ) 

Condition ( |2.3| ) translates into 

( [j k, a\ - ga^ f [JIT?-] ) x^ x'' + ( Ca - gai g"^ Cj ) 

r)V 8V 
+ ^d^a-9arr-d^)=^, a^l,...,m, (2.5) 
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In order to satisfy these equations it is sufficient to have 

Oaig"' [jk, r] = [jk, a] 

ga^rCr = Ca (2.6) 

dv _ dv 

These are the matching equations, see [ p^ . Following introduce variables 
relating the unknown mass matrix g to the original mass matrix g, 

K - 9a^ r ■ (2.7) 

Using A^, the matching equations take the form 

K b k, r] = [j k, a] 

K Cj = Ca (2.8) 

dV _ dV 

'dxi ~ 'dxP- 

Theorem 2.1. The following equations are equivalent to the matching equations 
A-equations: 

f) k — ^ r) 

— {g,,Xl)~[ka,2]\l-[kb,2]K^0, a,b=l,i.,m ^^.9) 



g-equations: 



v£ ^ , ^ I ^ 7? _ %j_ a=l,...,m 

" dx^ ^ dx^ ■ ^ dx^ ■ dx-' i,j^l,...,n 



F-equations: 



C-equations: 



■ dV dV 



Xi C, = Ca (2.12) 

For the proof see ||| , ||] . 

Remark 2.2. These equations always have a set of solutions of the form 
A^ = >^St g = -9 + 9", V ^ -V + V, Cj^ -C^ 

K K K 

with X ^ Q any constant, V° arbitrary function of the variables x'^ , £ ^ m + 1, . . . ,n, 
and g° any symmetric matrix valued function of the variables x^ such that g°^ = 0. 
We will call these solutions basic. 
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The A-equations are a system of ^ m(m+ 1) -n equations for n-m unknowns. It is 
not surprising that there are extra compatibiHty conditions. By viewing system (2.9) 
in the correct way, we are able to write down the compatibiHty conditions. Denote 



l^ab — 9ai ■ 



(2.13) 



Because the matrix gij is assumed to be non-degenerate, the matrix comprised of its 
m first rows has rank m. This imphes that out of m • n A's can be expressed as 
hnear combinations of z^'s, i.e.. 



A' 







- Vab- 

Substituting this in the A-equations, we obtain 

dk Vab - [a fc, a] h^^ vab ~ [b k, P] h^" v^a = [a k, p] A^ + [b k, p] Ag , 



(2.14) 



where index p varies over the remaining (n — m) indices. We wiU view system (2.14) 
of i m{m + 1) ■ n equations as a hnear algebraic system for the to (n — m) variables 



[^(/c,a,6)] K — F(kM,i 



(2.15) 



We know that this system has at least one solution by Rem ark |2.2| . Thus, the rank 
of the matrix A is at most m • n. In order for system ( [2.15 ) to have a solution, the 
vector F(^k^a,b) must be perpendicular to the kernel of the transposed matrix. A*, 



F _L ker A* . 



(2.16) 



Let the kernel of the matrix A* be generated by the vectors ■ The orthogonality 
condition (2.16) then takes the form of the following system of linear first order partial 
differential equations for Vab- 



^i'^'^^'^Hx) [dkVab - [ak,P]h^^Vdb - [bk, P]h^''vda] = 0. 



(2.17) 



Theorem 2.3. The general solution to the X-equations is given by any set of A^ 
solving the algebraic system {2. It), and A^ = h^"" Vab , where Vab is any solution to 
equations (2.11). 

In general, if to > 1, system ( ^.16 ) may be quite complicated and we do not have 
a satisfactory description of its solutions. 

3. Systems with one unactuated degree of freedom. If only one degree of 
freedom is unactuated, we do have a reasonable description of all solutions to system 
( 2.17 ). Assume, for simplicity, that gii{x) > 0. Then, after rescaling if necessary, 
we will have gii{x) = 1. More precisely, from the very beginning we could use, 
instead of (x^, x^, . . . , x"), the coordinates {z^, z^, . . . , z") which are related to x as 
follows: 

dz^ 



Vgiiix), 



In z coordinates the mass matrix is gij{z) ~ Qkiix) a nd, hence, gii{z) = 1. 

On the other hand, the structure of the equations of motion (2.1) does not change 
because of their tensorial form, and the condition ui = remains the same, again, 
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because ui = Uk ^p- = ui Thus, we assume that the coordinates are 

chosen appropriately, and gii{x) = 1. 

In the case of one unactuated degree of freedom one is solving for A^^. The A- 
equation reads 

^ = 2[kl,^]Xl, (3.1) 



where v — giiX\. Notice that [fcl, 1] = 0. View the equations (3T) as a system of 
linear algebraic equations for the variables A^', p — 2, . . . ,n. In order for this system 
of n equations in (n — 1) unknowns to have a solution, the vector 

must be perpendicular to the kernel of the matrix 

[11,2] ... [nl, 2]^ 



A* 



[1 1, n] ... [n 1, n] 



Let the kernel of A* be generated by the vectors £_r — (^r, ■ ■ • > O- The orthogo- 
nality condition for v translates into the system of equations 

X.(.)^C^(x)^ + ... + C(.)^ = 0. (3.2) 

The standard procedure to solve such a system of equations is to complete the system 
into an involutive system by adding equations [Xr, ^s]('^) = 0, [[Xr, Xg], XtKiy) = 
0, . . . , where [77* di, ('■' dj] = (?7'9i(C'^) — C* 9i{ri'^)) dk is the commutator of vector- 
fields. Recall that a system of equations 

Yiiiy) = 0,...,Yk{i^) - 

is involutive if [Yp, Yg] = fpg{x)Yr. 

Thus we have proved the following result. 

Theorem 3.1. With one unactuated degree of freedom there is a coordinate 
system such that the v-equations, (2.11), become a homogeneous linear system of 
equations for one unknown function. This system, (3.i ), can he completed into an 
involutive system. 

Remark 3.2. Note that we want to preserve the relationship (2.1), i.e., 

= gii - Xidji = 0. 
Using only the g- equations, one computes 

^ — + + ^gn + X"^ — - A^' ^ 

^ dx^ dx^ dx^ ^ dx^ dx^ 

i dEi dX{ „ d , g . g f dgei , dgu dgnX 



dx^ dx^ dx^ \ dx^ dx^ dx^ 
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Now invoke the X-equation to obtain 



0. 



We see that equality (2.1) holds locally provided it holds on a hypersurface transverse 
to X\. See also ^ Proposition 1.4] 

Given any non-zero solution of the A-equations, there is a local coordinate system 



y" such that 



(3.3) 



Let G and G represent g and g in the y-coordinates. The equation (2/7) then reads 

(3.4) 



dx"- 



The ^-equations read 



dGij dgki dx'^ dx^ 
dy^ dx^ dy^ dy^ ' 



(3.5) 



and the 1^-equations become 

dV _ dV 
dy^ dx^ 

It is easy to see that the following result holds. 

Theorem 3.3. Given any non-zero solution to the X-equation, there is a unique 
solution to the g- and V -equations with initial data prescribed at y^ = 0. 

Remark 3.4. Note that equation(3.4) gives directly 



Gki — 



ki 



d]f_ 

dx^ ' 



and so one only needs to solve ( jj*. 4 i for n{n — l)/2 quantities Gij, 2 < i, j < n. 

4. Example 1: Inverted pendulum in a vertical plane. As the first exam- 
ple we consider the inverted pendulum restricted to a vertical plane with horizontal 
and vertical actuation of the base, see Figure 1. 



O 




Figure 1 
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After rescaling units, the mass matrix and potential energy are given by 



1 —a cos(a;"'^) —a sin(x^)^ 

g = I —a cos(x^) 1 

-a sin(a;^) 1 

V = bx^ + cos(x^) 

Since only x^ is unactuated, we will simplify notation and use A' to denote A^^. The 
A-equations (2.9) are 

div = 2a sin(a;"'^) A^ — 2a cos(a;"'^) A^, = = 

with ^ = A^ — a cos(a;^) A^ — a sin(x^) A"^. It is not difficult to see that the general 
solution to these equations is 

Ai = v{x^) + - col{x'')dMx^) + 

2 sm(a;^) 



o ^VtT^i^(^') + coi{x^)\'{x\x\x'') 
la sm(x^j 



A — A ^ X , X ) , 

where v{x^), X'^ {x^ , x"^ , x'^) are arbitrary. In order to obtain a managable explicit 
solution to the matching equations, we will choose 



i^{x^) — a fiQ sin^(x"'^) + ao — a no, 
with free parameters (Tq and /ig. Then 

A^ = CTo , A^ = /io cos(x"'^). 

The coordinates 

1 1 1 2 2 -/In S 

y = — X , y ^ X - fio sm(x ) , y ^ 







satisfy (3.3). Following Remark 3.4, we need to solve the g-equations only for 322, 
g23, and 333. These equations are 

d ^ d ^ d ^ 

■^—[922 = TT^ 923 = TT^ 933 = 0- 

oy^ oy^ oy^ 

Clearly, 

?22 = 922{y^,y^) = 922{x'^ - Mo sin(a;^), x^) 
923 = 923{x'^ - Mo sin(a;^), x^) 
933 = 933{x^ - Mo sin(x^), x^) 
From g = 'g X, we obtain the rest of 'gtj as: 

I 



511 



+ ^ cos^x') + 4 cos'ix') 522 



Co 







(J, 







^ a MO ^ 1\ - 

ffi2 = cos(a; ) 522 

- a . / In Mo / i\ - 
513 = sm(a; ) cos(a; ) 523 

fQ CTo 
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The ^-equation yields 



The C-equation reads 
C = -oqR{x) 



V 



cos(xi) + w{y^,y^). 



0. One solution is 




ii^ COS(X^) 
1 
1 



CTo ^ ' 

1 

1 







I- 









The resulting control law can be obtained explicitly from equation (2.4). The expres- 
sion is too long to be included in this paper. 

Proposition 4.1. If the functions g22{y'^,y^), hsiv^.V^), dzziv'^.V^), w{y^,y^), 
and R(x), and the parameters and ctq are chosen so that 

?22(0) > 0, 323(0) = 0, 533(0) = 1, 

a^2w(0) > 0, dy2dy3w{0) = 0, a^3w(0) > 0, R{0) > 0, 
tJo < 0, 322(0) /io + a^iQ + ao > , 322(0) (a^o - fJo) + < , 

then X = X = is a locally asymptoitcally stable equilibrium of the closed loop 
system. 

5. Example 2: Inverted pendulum cart on a seesaw . In the previous 
example the kernel of the matrix A* , ( 2.15| ), was two-dimensional. Generically, for 
systems with one unactuated degree of freedom the dimension of the kernel will be 
1. The following example illustrates this situation. The inverted pendulum cart on a 
seesaw is shown in Figure 2. There are several interesting ways to actuate this system. 
We will consider the case with actuated cart and pendulum, and unactuated seesaw. 




Figure 2 

The rescaled mass matrix and potential energy of the system are given by 



a x^ sin(a;^ 







b + (a;3)2 

ax"^ sin(a;^ — x^) 1 —a cos(x-'^ — x^) 

—a cos(a;^ — x'^) 1 



and 



V = x^ sm{x'^) + a cos(a;^) 
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The theory in section 3 was presented with special coordinates so that gn — 1. 
However, in practice this is not necessary. 

As before, we write A* for A^^ and i' for gij X' . The A-equations are 

div = 2a a;^ cos(a;^ - x^) - 2x^\^ 

831^ ^ 2a sin(x^ - x'^) + 2x^ A^ 

Hence, v — i'{x'^,x^). Plug in A^ = {i' — gi2 X^ — gi^ X'^) / gn and solve for A^ and 
A^: 

Ai = 4(2^' - x'd^iy) 

^^- 2a5sin(l^-.^) (-^-^-+(^+(-^)^)^-- 



^^= 25.3sin|.^-.^) ('^(-^)^ cos(.^-.^). 

+ xH,H^r)cosix^-^^)d.^-^sH^^-x^)9.. 

Notice, that A^ and A^ blow up as x approaches unless 1/ = K{b+ (x^)^). 
Since g = g X, one must have detg ^ as x — > 0, i.e., g degenerates at x = 0. This 
means that H{x,x) = ^gijx'^x^ + V cannot serve as a Lyapunov function unless 
v = K{b + {x^)^). This v corresponds exactly to the basic solutions of the matching 
equations from Remark |]^. This illustrates the following general principle. 

Remark 5.1. // (xq, 0) is the desired equilibrium of a system and 

rankA*{xQ) < limsup ranfcj4*(a;) , (5-1) 

X — >xo 

then only basic solutions of the matching equations should be tested to produce a 



stabilizing control law from (2.4) 



6. Example 3: inverted pendulum cart on a roller coaster. Consider a 
cart with inverted pendulum on a roller coaster. Special cases of this mechanical 
system include the inverted pendulum on a rotor arm, the inverted pendulum on a 
verticle disk, and the inverted pendulum cart on an incline. By assuming that the 
size of the base of the cart is relatively small we may neglect the inertia of the base of 
the cart. It is therefore sufficient to model the cart with one point mass for the base, 
and one point mass a fixed distance away for the pendulum. The pendulum joint will 
be unactuated. 

The configuration of the system may be described by a position and an angle. As- 
sume that the shape of the roller coaster is given as a curve x{s) in parametrized 
by arc length, s, from a fixed point. Assume that the pendulum is always in the 
plane containing the tangent vector, t(s), and the vertical direction, 63. Let (f> be 
the angle between the pendulum and 63. By rescaling mass, length, and time, we will 
write 



9 
V 




b sinfa 



0) (l + fc(.)^sin^^(^^); 
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where a and b are positive parameters, < 6 < 1, and x is the vertical component 
of x(s) . The (miit) tangent vector to the curve is t(s) — , where ' stands 

for the derivative with respect to s. The curvature of the curve is k{s) = 
Denote by n(s) the principal normal to the curve. Recall that t'(s) = k{s)n{s). 
In the formula above is the vertical component of the principal normal, and a is 
the angle between r and the vertical direction. Index 1 corresponds to 0, index 2 
corresponds to s. The unactuated degree of freedom corresponds to the (j) variable. 



The A-equations (3J^) then read as follows: 

dii^ = —2b cos(q! — 0) A^ 



d2iy = k{sf sin(2 0) ( -) A? . 

sm a 



(6.1) 



The orthogonality equation (|3.2|) then, obviously, is 



k{sf sin(2 0) ( '^" ;^^ ) |5 + 2 6 cos(a - </>) |^ = . (6.2) 



sm a 



It is not clear if all solutions to this equation can be written explicitly for a general 
curve. We consider here two particular cases when this is possible. The first case 
is when sin^ a — rv^. This occurs exactly when the roller coaster lies in one vertical 
plane. The second case is when a{s) is constant. This occurs when the track is 
constantly inclined. 

6.1. Case 1: sin^ a — n\. Note that this case includes the interesting examples 
of an inverted pendulum on a vertical disk and an inverted pendulum cart on an 
incline. 



As is readily seen from ( |6.2D , the general solution of (6.2) in this case is = v{4>)^ 
an arbitrary function. Then 

2 _ 1 dv 

^1 — ^: 



and 



2 b cos(a 



1 djy 
Xl = v{(t>) + - tan(Q! - 0) — . 



This is a general solution of the A-equation. From here one must solve the g- and 
1^-equations. For special choices of a{s) and/or v{(j)) these equations have explicit 
closed form solutions. 

6.2. Case 2: a{s) = ao. Examples with a{s) constant include an inverted 
pendulum cart traveling on any path in a horizontal plane, a cart on a vertically 
oriented circular helix, or any constantly inclined track. 

Since ^ = — fc(s) , if a{s) — ao, then we have fc(s)ri3(s) = 0. To solve 
equation (^|^), we introduce new coordinates 

= f3{s) 

= /3(s) + cos(ao) In I CSC0 + cot0| — sin(Q;o) In I sec0 + tan(/)|. 



where 



lo bsm (ao) 
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Equation (6.2) is equivalent to the fact that v is an arbitrary function of . Thus, 
from (|6.l|) and the definition of v we find 



— v\z ) . i^.s — -r^ ; — 
sni(2(pj az^ 



1 di^ 
h sin(20) 



6.3. End of the roller coaster example. From the computations in case 1 
and case 2 one can see that the general solution to the matching equations for the 
cart on a roller coaster will be fairly complicated. However, we can show that any 
linear control law is the first order germ (linearization) of some matching control law. 
In fact, the only requirement for this is that there is no rank drop at the equilibrium. 



I.e., 



rankA*(a;o) = lim sup rank A* (x) 

X — >X^) 



(6.3) 



We assume that the dissipative term at the equilibrium satisfies the following natural 
assumptions: 



Q(xo,0)=0, ^ 



C, = 0. 



(2:0,0) 



Lemma 6.1. // condition \6.3^ is satisfied for a two degree of freedom system, 
then the first order germs of matching control laws at (xq, 0) exhaust all linear control 
laws for which the closed loop system has an equilibrium at (xo,0) . Proof. Given a 
linear control input 



lin 



Vj + Qij (x* - Xq) + bij X^ 



with v!^}'^ = 0, we will find a matching control law with the same germ. From the 
general expression (2^) for the matching control law, we see that the first order germ 
is 



f ™' = {V, - g,i g'"- Vr) + ( V,r - g,i % ){x''-xl)+ ( C,, - g,t C„ ) x' 
where 



dV 

dx^ 



dxWx^ 



Xq 



dO, 

dx^ 



ixo,Q} 



and Vj , Vjr , and Cji are defined similarly. Equating like terms gives 

V£] = da 5"' {Vr] - arj) , C'ij = gu g^'' {Crj - , = gu 5*'' (K- ~ Vr) ^ . 

(6.4) 



One can see that Vj, Vjr, and Cji are specified once gij{xo) is known. Moreover, 
there exists a non-degenerate symmetric gij{xo) such that the resulting Vij will be 
symmetric, see [Q, Lemma 1]. To conclude the argument, we now show that any non- 
degenerate symmetric ^ij{xo) arises as a zero order germ of a solution to g-equation. 
Also, any Vi , Vij satisfying Vi = 'gu g^^ {Vr — Vr) and Vij = gu g^^ (Vrj — arj) 
arises as a solution to the T^-equation. 
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Given a non-degenerate ^ijix^), define the non-degenerate A^(a::o) = gik{xo)d'^'' i^o)- 
Set 1^0 = gii{xo) X\{xo) + guixo) XI{xq). The A-equations in this case are 



div - 2 [1 1, 2] — 1/ 
511 



d2iy - 2 [1 2, 1] 



2 [12, 11 — 1/ 
511 



2 [1 1, 2] 
2 [21, 2] Xj 



(6.5) 



By the rank condition, we know that the rank of A* in the neighborhood of xq is either 
identicaUy or identically 1. If this rank is 0, then can be any constant times gn. 
We simply choose i^ix) = {i^o / gii{xo)) gii{x) ■ Any solution to the algebraic equation 
^{x) = gii{x) Xl{x) + gi2{x) Af (a;) is a solution to the A-equation. If the rank of A* 
is 1, then the orthogonality condition, ( ^.16 ), is 



[2 1, 2] diiy - [1 1, 2] d2iy + 2 ([1 1, 2] [1 2, 1] - [2 1, 2] [1 1, 2] ) — ly = . (6.6) 

511 

At X = xo, either di or 82 is not parallel to the vector [21, 2]9i — [11, 2] ^2 . 
Assume it is di . Then initial conditions to equation ( ^.6| ) can be specified along the 
line x^ — ^o- particular, we can choose the initial values so that 



aii/-2[ll, 2] — ly 
511 , 



2 [1 1, 2] 



A?(a;o), fixo) 



(6.7) 



X—Xq 



The second equation in (|6.5[), 



d2V - 2 [1 2, 1] — V 
511 , 



= 2 [2 1,2] 



X\{xo) , 



will be satisfied automatically since the rank of A* is 1. Let ^{x) b e a solution of 
(3.6) with initial condition, i/(a;^,0) , satisfying (6/7). Now, one solves (|6.5|) for Xl{x) 
and then i'{x) = gii{x) X\{x) + gi2{x) Xi{x) for A}(a::). 

Now that the A-equations are solved, we turn to ^-equations. These equations 
take the form 

■^g + Rg = S, 

where ~ X\di + Xf 82 ■ The initial conditions can be set on any line transverse 
to in particular, along the line X\{xq) (x^ — xj) + A^(a;o) {x^ — Xq) = 0. On 

t his l ine set A|(x) = A2(xo), and 5(2;) = g{x) ■ (A(xo))~^. The solution to equation 
(2.10) with this initial data then has the desired value ai x — xq. 

It remains to show that the F-equation has a solution such that — and 



Xiixa) Vij = Vrj - a. 



(6.8) 



where Af(xo) = grk{xo)d''^{xo)- Since Af(a;o) is non-degenerate, either Xl{xo) ^ 
or Af (xo) 7^ 0. Consider the case with A^(xo) ^ 0. In this case the line x^ = x"^ is 
non-characteristic for the ^-equation 



X\8iV + Xld2V ^ diV . 



(6.9) 
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Pick the initial value, V\^2-^2 , so that 



Vi = 0, Vn = 



and solve equation ([6.91 ). Since x = Xq is an equihbrium, Vi = and (V2 — "^2) = 0. 
From the differential equation, (|6^), we see that V2 = 0. Differentiating equation 
( |6.9D with respect to x^ and x^, we see that Wij = Vij satisfies 



By construction, 



A}(a:o)W^i2 + Xlixo)W22 ^ V12 ■ 



A2(a;o)W^ii + Xl{xo)Wi2 = F12 - ai2 



(6.10) 



(6.11) 



Notice that ( |6.4[ ) implies that Wgj — 'gug^'^ {Vrj — a^j) also satisfy equations ( 6.1Cl| ) 
and (6.11). Since the solution to the algebraic system (3.1C), (6.11) is unique, we 
conclude that equation (6.8) is valid, as required. 

7. Example 4: A double pendulum on a wheel. Our next example is the 
system with two unactuated degrees of freedom depicted in Figure 3. Only joint A is 
actuated. 

After rescaling, the entries gij of the mass matrix are 

gij — rriij cos(a;* — x-') 

and the potential energy is 

V — ai cos(x^) + a2 cos(a;^) + as cos(x^) . 
The parameters rriij = rriji and aj are positive. 




Figure 3 
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D. AUCKLY AND L. KAPITANSKI 



There are six unknown A^. Define Vab — dai '^b- Note that we must have 1^12 = i'2i- 
The computations in this section were performed using Maple. The A-equations ( |2.9| ) 
are 

di I'll = — 2 mi2 sin(2:i — X2) AJ — 2 mi3 sin(a;i — X2) Aj 

^2 i^ll = 

dsi^ii = 



di V22 = 

^2 i^22 = + 2 mi2 sin(a;i - X2) A2 - 2 77123 sin(a;2 - 2:3) \\ 

dzV22 = 



di yi2 = — 'mi2 sin[xi - X2) A2 - mi3 sin(a;i — 3:3) Aj 
^2 i^i2 = + '7112 sin(a;i - X2) A} - 77123 sin(a;2 — xs) A^ 
d-iVi2 = 

The second step is to express Aj, Aj, A^, and Aj, in terms of i^n, 2^12, and 1^22; 
After substitution into the above equations we obtain 

di vii = Dl 1 1^11 + Dl 1 7^12 + B\ 1 \\ 
82 vii = 
937^11 = 
dl V22 — 

d2V22 = £'3,2'^i2 + Dl 2V22 + Bl ^Xl (7.1) 

^3 7^22 = 

dl V12 = Dl 1 7^12 + Dl 1 7^22 + ^2^1 A2 
82 V12 = ^11 + ■Di_2 ^^12 + -82,2 ^1 

^3 i^l2 = 

Here the D^ ^ and are expUcit expressions involving x. In Section 2 we de- 

scribed a general procedure to obtain compatibility conditions for this system. In 
this particular case, however, we use a different tactic: we compute and compare the 
mixed derivatives of Vab- The first set of equations we obtain is 

dzdi vi2 = Kii 93Ai + A'i2 A^ = 

53^2 V22 = ^^21 d^\l + i^22 = 

Direct computation shows that AGt[Kij) ^ 0. Hence, Af = 0. Similarly, 

d^di vii = Lii dsXl + L12 Xl = 
9382 V12 = ^21 ds,\\ + L22 Aj = 

and dei{Lij) ^ 0. Hence, \\ — 0. 

Next, we substitute Af = A2 = into (7.1) and solve for j/n, 7/12, 7^22- This 
gives 

, ^22 mi2 2 In 

7^11 ~ const, 7^22 = Vll , V12 = COS(X — X ) Vll 

mil mil 
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Returning to A-equations we see that 

Ai = A9 - z/11 , An = Ai = . 

mil 

Our computation shows that the only solutions of the matching equations are the 



basic solutions defined in Remark 2.2 
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